Observation of Anti-correlation of Classical Chaotic Light 
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We wish to report an experimental observation of anti-correlation from first-order incoherent 
classical chaotic light. We explain why the classical statistical theory does not apply and provide a 
quantum interpretation. In quantum theory, either correlation or anti-correlation is a two-photon 
interference phenomenon, which involves the superposition of two-photon amplitudes, a nonclas- 
sical entity corresponding to different yet indistinguishable alternative ways of producing a joint- 
photodetection event. 

PACS numbers: PACS Number: 03.65.Bz, 42.50.Dv 



In 1956, Hanbury Brown and Twiss (HBT) discovered 
a nontrivial intensity correlation in thermal light [l[ . Fig- 
ureQ]schematically illustrates a modern HBT interferom- 
eter or the so called "intensity interferometer" . In a tem- 
poral HBT interferometer, the temporal, randomly dis- 
tributed, chaotic thermal light has a twice greater chance 
of being measured within its coherence time by the joint- 
detection of two individual photodetectors. In a spatial 
HBT interferometer, the spatial, randomly distributed, 
chaotic thermal light exhibits a twice greater chance of 
being captured within a small transverse area that equals 
the spatial coherence of the thermal radiation by two 
point photodetectors . It was recently found that for a 
large angular sized chaotic source the spatial correlation 
is effectively within a physical "point" . The point-to- 
point near-field spatial correlation of chaotic light has 
been utilized for reproducing nonlocal ghost images in a 
lensless configuration 0]. 
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FIG. 1: Schematic of a modern HBT interferometer which 
measures both temporal and spatial correlation of light by 
scanning the optical fiber tips longitudinally or transversely. 



What is the physical cause of this peculiar behavior 
of chaotic thermal light? In the classical point of view, 
the HBT phenomenon is a statistical correlation of in- 
tensity fluctuations of thermal radiation. In general, no 
matter how complicated the optical setup is, classical the- 
ory considers the joint-detection between two individual, 
point-like photodetectors, D A and Db, as a measure of 
the statistical correlation of two intensities at space-time 



coordinates (r A ,t A ) and (r^,^) 

T^(r A ,t A ;r B ,t B ) = (I(r A , t A ) I(r B , t B )) (1) 
= (I(r A ,t A )}(I(r B ,t B )) + (M(r A ,t A )AI(r B ,t B )). 

The HBT observation comes from the second term of 
Eq. |T]) , which indicates correlated intensity fluctuations 
at a distance. Why do the two distant intensities fluc- 
tuate in such a peculiar manner? One historical answer 
is "photon bunching", i.e., a thermal light source has a 
higher chance of emitting photons in pairs. Another in- 
volves the coherence of the electromagnetic fields 

T^(r A ,t A ;r B ,t B ) = (I(r A , t A ) I(r B , t B )) 
= (E*(r A ,t A )E(r A ,t A )E*(r B ,t B )E(r B ,t B )) 



_ r (i) r (i) 
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where 



defined as the self-coherence func- 



J = A,B, 
tion, or self-correlation function of the field; Y 



(i) _ 

AB 



^*ba = (E*{r Al t A )E(r Bl t B )} is defined as the mutual- 
coherence, or mutual-correlation function of the field. It 
has been well accepted that an intensity-interferometer 
measures the mutual-correlation T^ B of the input field, 
no matter how complicated the optical setup. 

4a= iIb 




FIG. 2: Schematic setup of the experiment. 

Now we consider a slightly different experimental setup 
in Fig. [21 Instead of directly measuring the intensity 
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correlation of (I{r A , t A ) I{r B} ta)), this setup measures 
(/(ri,ti) I(r 2 , t 2 )) at the output ports of a 50/50 optical 
fiber beamsplitter. If the two input fiber tips A and B 
are placed within the longitudinal coherence time and the 
transverse coherence area of the thermal field, this setup 
is equivalent to a Mach-Zehnder interferometer. D\ and 
D 2 will each measure first-order interference as a function 
of the optical delay S when scanning the fiber tip A along 
its longitudinal axis. Consequently, the joint-detection of 
D\ and D% outputs an interference pattern that is fac- 
torizable into two first-order interferences. However, this 
is not our desired experimental condition. We decided to 
move the fiber tip A outside the transverse coherence area 
to force = 0. Under this condition, there would be 
no observable first-order interference in a timely averaged 
measurement and in a instantaneous "single-exposure" 
observation . Consequently, the standard classical sta- 
tistical correlation theory of intensity fluctuations gives 
the following second-order correlation T^ 1 that is inde- 
pendent of optical delay <5, 



(3) 



rP)(ri ) ti;r 2s *2) = h) I(r 2 , t 2 )> 

= (\E(r A ,t 1A ) + E(r B ,t 1B )\ 2 \E(r A ,t 2A ) - E(r B ,t 2B )\ 2 ) 

where (r^,*^ = tj - r A j/c) and {r B ,t 3B = tj - r Bj /c), 
j = 1, 2, with r A j/c (r B j/c) the optical delay from point 
A (B) to Dj, are earlier space-time coordinates at points 
A and B. In the Appendix A we show in detail why the 
classical statistical intensity correlation r( 2 )(ri, t\ \ r 2 , t 2 ) 
is 5 independent. In fact, the physics is rather simple; 
mathematically Eq. ((3]) has sixteen terms to begin with. 
However, all ^ AB l ab terms vanish due to the mutual in- 
coherence between fields E(r A ,t A ) and E(r B ,t B ), leav- 
ing only these nonzero terms which have no contribution 
from the 5-dependent first-order mutual coherence. 

The measurement produces quite a surprise. A "unex- 
pected" anti-correlation "dip" was observed in the joint- 
detection counting rate of D\ and D 2 as a function of 
the optical delay S. Figure [3] reports two typical mea- 
sured anti-correlation functions with different spectrum 
bandwidths of the chaotic field. 

The experimental detail is described as follows. 

(1) The source: the light source is a standard pseudo- 
thermal source that was developed in the 1960's and used 
widely in HBT correlation measurements Q ■ The source 
consists of a CW mode-locked laser beam with ^200 fem- 
tosecond pulses at a 78 MHz repetition rate and a fast 
rotating diffusing ground glass. The linearly polarized 
laser beam is enlarged transversely onto the ground glass 
with a diameter of 4.5mm. The enlarged laser radiation 
is scattered and diffused by the rotating ground glass to 
simulate a near-field, chaotic-thermal radiation source: a 
large number of independent point sub-sources with in- 
dependent, random, relative phases. 

(2) The interferometer: A 50/50 beam splitter (BS1) 
is used to split the chaotic light into transmitted and re- 
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FIG. 3: Two typical observed anti-correlation functions with 
different coherence time of the chaotic field. r c ~ 345/s for 
(I), t c ~ 541/s for (II). The coherence time is determined by 
the bandwidth of the spectral filters (IF). 



fleeted radiations which are then coupled into two iden- 
tical polarization-controlled single-mode fibers A and B 
respectively. The fiber tips are located ~ 200mm from 
the ground glass, i.e. d A — d B ~ 200mm. At this dis- 
tance, the angular size Ad is ~ 22.5 mili-radian (1.29°) 
with respect to each input fiber tip, which satisfies the 
Fresnel near-field condition. Two identical narrow-band 
spectral filters (IF) are placed in front of the two fiber tips 
A and B. The transverse and longitudinal coordinates of 
the input fiber tips are both scannable by step-motors. 
The output ends of the two fibers can be directly coupled 
into the photon counting detectors D\ and D 2l respec- 
tively, for near-field HBT correlation measurements, or 
coupled into the two input ports of a 50/50 single- mode 
optical fiber beamsplitter(BS2) for the anti-correlation 
measurement. 

(3) The measurement: two steps of measurements 
were made. The purpose of step one is to confirm the 
light source produces chaotic-thermal field. We mea- 
sured the HBT temporal and spatial correlation by scan- 
ning the input fiber tips longitudinally and transversely. 
In this measurement the output ends of the fibers are 
coupled into D A and D B directly as shown in Fig. Q] 
Chaotic radiation can easily be distinguished from a laser 
beam by examining its second-order coherence function 
G' 2 ' (r A , t A ; rs, t B ), which is characterized experimen- 
tally by the coincidence counting rate that counts the 
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joint photo-detection events at space-time points (r^, tjCj 
and (i\b,£,b). Figure [4] reports two measured second- 




FIG. 4: Measurement of G^ 2 '(xa — xb) at no rotation(0 rpm) 
and at 1000 rpm. Here, xa and xb are the x-component of qa 
and pB, and correspondingly y components are kept tja — ys- 

order correlations at zero and at 1000 revolutions per 
minute (rpm) of the rotating ground glass. This mea- 
surement guarantees a typical HBT correlation of chaotic 
light at rotation speeds greater than 1000 rpm of the 
ground glass, indicating the chaotic nature of the light 
source. In this measurement, we have also located ex- 
perimentally the longitudinal and transverse coordinates 
of the fiber tips A and B for achieving the maximum 
coincidence counting rate, corresponding to the maxi- 
mum second-order correlation. In step two, we couple the 
50/50 fiber beamsplitter (BS2) into the setup as shown 
in Fig. [21 This measurement was done in two steps. We 
first measured the first-order interference at pa = Pb 
by scanning the input fiber tip A longitudinally in the 
neighborhood of cLa ~ 200mm. There is no surprise to 
have first-order interference in the counting rates of D\ 
and D 2 respectively. When choosing pa = Pb-, the two 
input fiber tips are coupled within the spatial coherence 
area of the radiation field; we have effectively built a 
Mach-Zehnder interferometer. We then move the input 
fiber tip A transversely from pa — Pb to \pa — Pb\^> l c , 
where l c is the transverse coherence length of the thermal 
field. Then we scan the input fiber tip A again longitu- 
dinally in the neighborhood of <Ia ~ 200mm. The opti- 
cal delay between the plane z = d,A ~ 200mm and the 
scanning input fiber tip A is labeled as 6 in Fig. [2] We 
have thus achieved the expected experimental condition 
of r^L = 0. There is no surprise we lose any first-order 
interference in this experimental condition. However, it 
is indeed a surprise that in the joint-detection of D\ and 
D 2 an anti-correlation is observed as a function of the 
optical delay 5 that is reported in Fig. [3l In these mea- 
surements, \pA — Pb \ ~ 40Z C . 

Quantum theory gives a reasonable explanation to this 
surprising observation. In quantum theory, the second- 
order correlation function represents the probability of 



observing a joint-photodetection event at space-time co- 
ordinates (ri,ii) and (r 2 ,i 2 ) 5]. If more than one, dif- 
ferent yet indistinguishable, alternative ways of trigger- 
ing a joint-photodetection event exist, these probability 
amplitudes must be linearly superposed resulting in an 
interference. The HBT correlation of chaotic light is the 
result of two-photon interference, which involves the su- 
perposition of two-photon amplitudes, 

G^(v 1 ,t 1 ;r 2 ,t 2 ) (4) 
= tr [ p (n , t x )#-> (r 2 , ta)#+) (r 2 , t 2 )E i+) (r x , h ) ] 

= ^P a p\-^=[^{Y a i,ti;vp2,t 2 ) +^(r / 3i,ti;r a2 ,*2)] | 2 , 

a,/3 

where E^(rj,tj), j = 1,2, is the positive (+) or nega- 
tive (— ) field operator at coordinate (rj,tj). In Eq. ([4]), 
we have treated the chaotic radiation in a mixed state 

P ~ |Q)(0| +£)P Q |* Q )<* Q | +2^l*a>l*/»X*/»K*a| + 
a a,/3 

which represents an ensemble of sub-radiations, such as 
trillions of photons created from a large number of inde- 
pendent and randomly radiated sub-sources. In the third 
term of p, which contributes to the joint-photodetection, 
Pa/3 represents the probability for the ath and the /3th 
sub-radiations to be in the states 1^)1^^). We have also 
defined an effective wavefunction G^l = \^ a p\ 2 with 

* Q/ 3 = -j=[s^(r al ,h;ri3 2 ,t 2 ) + d{rpi,ti ;r a2 ,t 2 )] (5) 

where s/(r al , h; r p2 , t 2 ) = (0\e[ +) |* a ) (0|^ +) |^) [or 
^(r^i,ti;r a2 ,t 2 ) - (0|^ +) |^)(0|^ +) |* Q )] represents 
the alternative in which the ath [or /3th] photon and the 
/3th [or ath] photon are annihilated at (ri, tx) and (r 2 ,t 2 ), 
respectively. The symmitrized effective wavefunction in 
Eq. ([5]) plays the same role as that of the symmitrized 
wavefunction of identical particles @ . Obviously, it is the 
superposition of srf(r al , t x ; r 02 , t 2 ) and ja/(rpi,ti; r a2 ,t 2 ) 
causing the nontrivial HBT correlation. 

In the view of quantum theory, anti-correlation is ob- 
servable from classical chaotic light, if 

*a/3 = -^=[^(r Q i,ti;r^ 2 ,i 2 ) - sf(rpi, tf, r a2 , t 2 )] (6) 

is achievable. A number of experimental approaches may 
achieve the two-photon destructive interference condition 
of Eq. In this experiment, it is the beamsplitter 

(BS2) that breaks the symmetry of the effective wave- 
function in Eq. ([5]) and introduces the "— " sign in achiev- 
ing Eq. ([6]). In fact, both destructive ("— ") and construc- 
tive ("+") two-photon interferences have been achieved 
in the measurement of entangled two-photon states jjjj]. 
The anti-correlation "dip" has been widely used for iden- 
tifying nonclassical states since the mid 1980's [8(. 
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In the view of quantum theory, either the correlation 
or the anti-correlation of thermal radiation is the result 
of two-photon interference, which involves the nonlocal 
superposition of two-photon amplitudes Analogous 
to Dirac's statement that a photon interferes with itself, 
this interference is a jointly measured pair of indepen- 
dent photons interfering with the pair itself. Figure [5] 
schematically illustrates four alternatives for two inde- 
pendent photons to trigger a joint-detection event of D\ 
and D 2 . In (a) and (b) the measured pair comes from 
the same fiber tip, A or B. In (c) and (d) the measured 
pair comes from different fiber tips, one from A and the 
other from B. It is the superposition between amplitudes 
in (c) and (d) that produces the anti-correlation. 
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FIG. 5: There are four alternative ways for a measured pair 
of independent photons to trigger a joint-detection event of 
Di and D2. 

Now we give a formal quantum mechanical calcula- 
tion. The probability of observing a joint-detection event 
at space-time coordinates (ri,ii) and (r 2 ,t 2 ) is given in 
Eq. ([4]), in which the fields E(ri,ti) and E(r 2 ,t 2 ) are 
treated as the superposition of the early fields at the in- 
put planes z = d,A and z = oIb 

E(+\z 2 ,t 2 ) = -L[#+)( TA2 ) _ £(+)( TB2 )]. (7) 

where r A j = tj — rAj/c, j = 1, 2, the optical delay from 
the detector Dj to the input plane A, n the index of 
refraction of the fiber, and tbj is defined similarly. 

To simplify the calculation, we assume two groups of 
randomly distributed wavepackets are excited at the fiber 
tips A and B 

= / dwf{u)e^ A ot(«)|0) 
|tf B ) = / duf{u>)e^ B b\u)\Q), ( 8 ) 



where /(w) is the spectrum function which is mainly de- 
termined by the spectral filters (IF) in the experimen- 
tal setup, to a and t^B represent the times when the 
wavepackets pass the planes A and B, respectively. 



Substituting the field and the density operators into 
Eq. ([4}, it is straightforward to find that 



G {2) (z 1 ,t 1 ;z 2 ,t 2 ) 
+ k(TJi,Tl 2 )-^(rf 2 ,Tf 1 )] 



R , 



B2) 



(9) 



with 



^(tL.t&J - (0\E^(t a1 )\* a )(0\E^(t A2 )\V a ) 
) = <0| J B (+) (r B i)|* B )(0|£ ;(+) (T B2 )|* B ) 
(0\E(+\t a1 )\V a )(0\E(+\t B2 )\* b ) 



T l7 T B2) 

(r#i,rf 2 ) = (0\E( + Ht b1 )\*b)(0\E( + \t A2 )\*a), 



where, t\\ = {ti-t 0A )-r A x/c (rj^ = (ti-t 0B )-r B i/c) 
indicating the early time at points A (B) by propa- 
gating through the reflected (transmitted) optical path 
from source A (B) to the photo-detector D\\ and t^ 2 = 
(*2 -toA)- r A2 /c (t£ 2 = (t 2 -t 0B )- r B2 /c), indicating 
the early time at points A (B) by propagating through 
the transmitted (reflected) optical path from A (B) to 
the photo-detector D 2 ; to a and £ob are the initial times 
at A and B, respectively. These four effective wave func- 
tions correspond to the four alternatives shown in Fig|5] 
It is not difficult to see that the first two terms in Eq. (J9j) 
are 5 independent. The two-photon interference occurs in 
the third term. The cross two-photon interference term 
in the third term gives 



^*(ri 



l! 



for a Gaussian spectrum function of f(u>). The details are 
given in Appendix B. The coincidence coincident count- 
ing rate is therefore, 



R c ~ R J dt\dt 2 G ( - 2 \zi, t\\ z 2 , t 2 ) cx 1 - -1 



£2 /_2 



(10) 

Eq. pop has been verified by achieving different coher- 
ence time of the chaotic filed, which are determined by 
different bandwidth of the spectral filters (IF), as shown 
in Fig[3](I) and (II), except with lower contrasts [~ 28% 
for (I) and - 29% for (II)]. 

In conclusion, we have observed a nonclassical anti- 
correlation from classical chaotic light under the exper- 



imental condition of G 



(i) 

AB 



0. The classical statisti- 



cal correlation theory seems unable to explain this ex- 
perimental result. This observation is different from all 
historical measurement of the "dips" , either observed 
from nonclassical sources or from synchronized coher- 
ent sources. In the view of quantum mechanics, cither 
the anti-correlation "dip" or the correlation "peak" of 
thermal light are straightforward two-photon interference 
phenomena, involving the superposition of two-photon 
amplitudes. 
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APPENDIX A: CLASSICAL STATISTICAL 
INTENSITY CORRELATION 

The standard statistics of classical chaotic light gives 
the second order correlation in Eq. |3j). Eq. (J3j) can be 
explicitly calculated as 

T^(r u t i; r 2 ,t 2 ) = (/(r^ii)/^^)) 

= (\E(r A ,t 1A ) + E{r B ,t 1B )\ 2 \E(r A ,t 2A ) - E(r B ,t 2B )\ 2 ) 

= (\E(r A ,t 1A )\ 2 ) (\E(r A ,t 2A )\ 2 ) 

+ (E*(r A ,t 1A )E(r A ,t 2A )) (E(r A , t 1A )E*(r A , t 2A )) 

+ (\E(r B ,t 1B )\ 2 ) (\E(r B ,t 2B )\ 2 ) 

+ (E*{r B ,t 1B )E(r B ,t 2B ))(E(r B ,t 1B )E*{r B ,t 2B )) 

+ (\E(r A ,t 1A )\ 2 ) (\E{r B ,t 2B )\ 2 ) 

+ (\E(r B ,t 1B )\ 2 ) (\E(r A ,t 2A )\ 2 ) (11) 

- (E* (r A ,t 1A )E(r A ,t 2A )) (E(r B , t 1B )E* {r B ,t 2B )) 

- (E{r A ,t 1A )E*(r A ,t 2A )) (E*(r B , t 1B )E^r B , t 2B )) 

where, again, (r A ,tj A — tj — rj A /c) and {r B ,tj B — tj — 
rj B /c), j — 1,2, are the earlier space-time coordinates. 
Eq. (fTT|) can be written in the following form 

r (2)_ r (l) r (l) , r (l) p(l) , r (l) p(l) , r (l) r (l) 
L ~ L All L A22 "r" L All L A21 L Bll 1 B22 ' 1 B12 L B21 

, r (i) r (i) , r (i) F (i) _ r ( 1 ) pM _pW pW 

"r" 1 .411 L B22 "r" 1 A22 1 Bll L A12 L B21 L A21 L B12' 

(12) 

Here, all the terms T^r^ vanish due to mu- 
tual incoherence nature between the fields E(rA,tA) 
and E(r Bl t B ), leaving the non-zero terms T^r^, 

T BB r BB and T< AA r BB contributing to . The nega- 
tive signs in the last two terms are basically introduced by 
beamsplitter. It is not too difficult to find from Eq. (jTTJ) 
that r^ 2 - 1 is independent of the optical delay 5 for either 
CW or pulsed chaotic light. 

Case (I): CW chaotic light 

There is no doubt the first four terms are S inde- 
pendent. In these terms either the self-correlations, 
-^aii-^A22 an( i r^^r^^; or the cross-correlations, 
T^r^i and r^r^i are respectively associated with 
the same radiation A or B. The other four terms may 
contain 5 in their amplitude part or in their phase part. 
We do not need to worry about the amplitude part 
due to the stationary nature of the CW chaotic light. 



Let us examine the phase part: (1) the first two terms 

^lli^sL = {IauHIbv) and r^rj& = (/A22>(/Sll) 
contain the expectation of intensities only and thus are 

phase independent. (2)The second two terms r^j^r^i 
and r^iTs^ may contain relative phases of k(r A \ — r A2 ) 
and k(r B i — r B2 ), however, these phases are 5 indepen- 
dent. We may conclude in the case of CW chaotic light, 
I^ 2 ) is independent of the optical delay 6. 

Case (II): Pulsed chaotic light 

All the analysis are the same as above, except we 
do need to consider the delays between the amplitudes, 
which is 8 dependent. In the pulsed case, the fields are 
non-stationary, the last four terms have nonzero values 
only when their amplitudes have nonzero values simulta- 
neously in the joint-detection of D\ and D 2 . It is clear 
when (5 = 0, all four terms achieve their maximum values 
simultaneously due to the complete overlapping of their 
respective amplitudes. The overall contribution of the 
four terms, however, is null to due to the cancela- 
tion between the first two positive contributions and the 
last two negative contributions. When 6 increase or de- 
crease from zero, although the overlapping between the 
A amplitude and the B amplitude becomes smaller and 
smaller yielding smaller contributions to , the overall 
contribution keeps its zero vale for the same reason. We 
may conclude that in the case of pulsed chaotic light, T^ 2 ' 
is independent of the optical delay 6 too. 



APPENDIX B: QUANTUM THEORY 

In Appendix B we show a simple calculation of the ef- 
fective wave functions. The effective-wavefunctions cor- 
responding to the case in which the joint-detection event 
is produced by two photon coming from the same point 
is shown below as an example 

<*lTAiA) = (0|£ (+ W)I*a)(0|£(+>(ta 2 )I*a) 
oc J dujf{uj)e- iuJT ^ J dcj'f{uj')e~ tuJ ' T ^ 

= e-^^ifi^e-^^ {/(*/)} (13) 
where, 

T T {f{u)} = J duf{v)e- ivr 

is the fourier transform of the spectrum function f(v), 
and v = uj — luq (ujo is the central frequency) . It is easy 
to see these functions are S independent. 

The superposition of the two-photon amplitudes (c) 
and (d) corresponding to the interference term in Eq. ^ 
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is calculated to be 

HtI^-^t^t*)! 2 (14) 

= I^l{/M}^ 2 {/M}| 2 + l^{/M}^ 2 {/M} 

It is easy to find that the first two terms in Eq. (fl4|) 
contribute a constant to the coincidence counting rate 
R c . The nontrivial contributions come from the last two 
cross terms. Assuming a Gaussian spectrum, the cross 
terms is approximately to be 

for synchronized radiations A and £>, i.e., toA = £os, 
where r c is the coherence time of the measured field. It 
is interesting to see the cross term contribution is time 
independent. The coincidence counting rate is therefore 

R C ~R j dt 1 dt 2 G {2) (z 1 ,t 1 ;z 2 ,t 2 ) oc 1 - \^ !t '- 
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